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The evolution of spin network states in loop quantum gravity can be defined with 
respect to a time variable, given by the surfaces of constant value of an auxiliary scalar 
field. We regulate the Hamiltonian, generating such an evolution, and evaluate its 
action both on edges and on vertices of the spin network states. The analytical 
C . computations are carried out completely to yield a finite, diffeomorphism invariant 

result. We use techniques from the recoupling theory of colored graphs with trivalent 
vertices to evaluate the graphical part of the Hamiltonian action. We show that the 
action on edges is equivalent to a diffeomorphism transformation, while the action on 
' vertices adds new edges and re-routes the loops through the vertices. A remaining 

usresolved problem is to take the square root of the infinite-dimmensional matrix 
of the Hamiltonian constraint and to obtain the eigenspectrum of the "clock field" 
Hamiltonian. 
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I. INTRODUCTION 



o 

In canonical quantum gravity the notion of evolution requires careful definition because the trans- 
lations in time direction can be interpreted as diffeomorphism transformations Q. To be able to 
talk about evolution we can use "relational constructions" . Some physical fields can be introduced 
as a reference frame with respect to which the evolution can be defined Q], || , Q . In loop quantum 
gravity || Rovelli and Smolin || define a time variable by the surfaces of constant value of 
an auxiliary scalar field. By fixing a gauge in this construction the infinite number of Hamiltonian 
?H . constraints (one per space point) reduces to one constraint which can be interpreted as a Schrodinger 

equation and a Hamiltonian operator can be identified. We use the Hamiltonian obtained in this 
model to investigate the evolution of quantum gravitational states. In loop quantum gravity spin 
network states, JM, |§f|, ^|, furnish a complete basis of quantum kinematical states. 

There exists pQ ] , [ |ll| already a well established procedure for expressing different quantities from 
quantum gravity in terms of the loop variables frjjfl . Then operator versions of the gravitational 
quantities can be defined by replacing the loop variables with the corresponding operators. Thus it 
is relatively straightforward to introduce a loop version for the Hamiltonian operator. We use the 
result from [Q and 0] as a starting point for our calculations. Our purpose is to determine in detail 
the way the spin network states evolve under the Hamiltonian introduced in As we will see the 
result of the evolution can be split into two parts. First there is a multiplicative factor which is finite 
and diffeomorphism invariant. Second - the spin network states evolve topologically: The new state 
is a sum of terms, each term being based on the original spin network with an added extra edge of 
color 1. The added edge connects pairs of the original edges, meeting at a vertex. Also some change 
of coloring of the edges occurs such that the new graph is a spin network again. 

The content of the paper is as follows: In section 2. we define the Hamiltonian and regulate it 
to show that it has well defined action on the spin network states. We also introduce a modifica- 
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tion in the way the loop operators are defined, better suited for our calculations. Technically the 
action of the Hamiltonian operator can be split into analytical and graphical parts. The analytical 
part includes various pre-factors and integrals. The graphical part expresses the topological trans- 
formations occurring in the spin networks. In section 3. we compute the analytical action of the 
Hamiltonian separately on edges and on vertices. We show that to a great extent the action on 
edges is equivalent to diffeomorphism tranformation. Using some techniques from recoupling theory 
of knots and links with trivalent vertices we perform the graphical computation of the action of the 
Hamiltonian in section 4. The result from the graphical calculation tells us whether the diffeomor- 
phism class of the spin network or the coloring of certain edges change. We conclude with discussion 
of some open issues. 



II. THE HAMILTONIAN OF THE THEORY 

Because of the absence of external time with respect to which the evolution can be defined, we 
need some additional construction. As it has been shown in j| to define time we can use the physical 
degrees of freedom of an auxiliary field. We start by introducing a scalar field T(x). To serve as a 
clock, this field should be monotonically increasing everywhere on the space manifold sQ. Then we 
can use its 3-surfaces of constant value T(x) — const to represent the time with respect to which 
the evolution will be defined. The scalar "clock" field can be incorporated in the theory through the 
standard Klein-Gordon Lagrangian: 



£t = ^(x)y/^g(x)d^T(x)d u T(x). (1) 

In this expression [i plays the role of a coupling constant between the scalar field and the grav- 
itational field. Later on we will treat T(x) as time so, from dimensional analysis, the constant fx 
should have dimensions of energy density. The momentum conjugate to the field T(x) will be: 

clC 

* ( * ) = 0^ = M ^ ^ (2) 
Performing a Legendre transformation we get for the total Lagrangian: 

~ 2 

C = nd T - N^- - N^q 2 q ab d a Td b T - nN b (d b T) + C Gramty (3) 

The Euler-Lagrange equations we obtain from this Lagrangian are: 

S -^ = d a T-N-- N\d b T) = 0, (4) 
oir n 



and 



sc 



-d n + nd b (N f C d c T) + d b {N b n). (5) 



At this point we impose a gauge fixing, restricting the freedom of choosing the time coordinate. 
The gauge we use is d a T(x) = which, because of Eq. (0), implies that: 



d a T = N-. 



1 We use the standard notation - g^ v and q ab are respectively the 4-metric on M = S x R and the 3-metric 
on E. a,b, . . . are spatial and - internal indices; they all run from 1 to 3. 
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Thus the lapse function N(x) should satisfy the relation: 

N(x) = —— -, 6 

TT(X) 

where a(t) is an arbitrary function of (the coordinate) time. For the Hamiltonian constraint we get: 

C 0) = ^+ Cg (x), 

where 

Cg (x) = e ijk E ai E b ^F^ - Aq = C Emstem - Aq 

is the gravitational Hamiltonian constraint in terms of the Ashtekar variables and A is the cosmo- 
logical constant. According to the general prescription for gauge fixing in constraint systems we 
have to compute the Poisson bracket between the gauge and the existing constraints to check for 
secondary constraints. Thus we get: 

{C (x'),d a T(x)} = d a {^^-,T(x)} = d a [^-S 3 (x',x)}. 

The above expression vanishes when we smear the Hamiltonian constraint C (x 1 ) using N(x) from 
eq.Eq. (Q). Thus the only constraint which remains to be imposed on the wave functionals of the 
theory is the integral of the Hamiltonian constraint with the lapse function: 



Aff W c (.) = «(*/^ + V^&W) - V^&W) - o. 

Because we are going to impose this integral as a constraint operator, we can think of tt(x) as 



being equal to y — 2/i Cg (x) Thus the expression in the first parentheses can be replaced by ^/ ^tt(x) 
and the whole integral reduces to: 



a{t) J d 6 x \ jr(x) - y~2fiC G (x) J S 0. (7) 

Note that we assume that the Hamiltonian constraint Cg(x) satisfies the weak energy condition 
jyl, which in this case requires that Cg(x) < 0. 

In the process of quantization we promote this constraint into an operator equation. We require 
that in the loop representation the spin network states, depending also on the clock variable T, are 
annihilated by the constraint operator: 



(S,T\a(t) j ,f'.r | rr(.r) V -2//CT f> -(.r) | 0. (N) 
We interpret the integral: 



d 3 XTr(x) 

s 



as ill times a derivative with respect to the "clock" field and thus we arrive at the Schrodinger 
equation: 
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ih-^{S,T\ = {S,T\H (9) 



where the Hamiltonian H corresponds to the classical observable fjfwe have restored the constants) 



H = -^E= / <Pxy/-C(x)+Aq{x). (10) 



s 



where 

- C(x) + Aq(x) = -e ijk ^ b E<"E b i - ±Ae abc E<" E b > E ck } . (11) 



We will omit the overall factor of (^fjl / V 'AttG) in what follows, as it is not important for our 
considerations. 

Because of the product of distributional quantities in the Hamiltonian we have to regulate the 
above expression. Using an arbitrary, e.g. Euclidian, background metric h ab (x) we divide the space 
manifold £ into cubes of size L, labeled i?/. In the calculations we will let L to go to zero as the 
number of cubes goes to infinity. Then we can write Eq. © MS 



H = lim f [ d^xVh{x-)J - C{X l+ Aqix) = lim 



j3 „./r^./- c W + A «W _ sr^.l-^ + ^il, (12) 



h{X!) 



where h(x) is the determinant of the background metric, and xi in the last expression is a point in 
the cube Ri. Now putting L 3 inside the square root and going back to integral form we get: 



/ r d^r i 

H = I^E^ 3 J Ri y^y (-^(*) + M*)) - limE\/-^ + A ^ ( 13 ) 

All of the above manipulations are correct in the limit L — > 0. The last step in Eq. (|l3j) is a definition 
of C/. Vf is the classical volume of the J-th cube. 

Let us now focus on the first term of the expression under the square root in equation Eq. Jl3|). 
We introduce a regulating point-splitting function fs{x,y) which is a density of weight one with 
respect to its first argument. This function satisfies the requirement that for any smooth function 
<t>{x): 

lim / d*xct>{x)f s {x,y) = <j>{y). (14) 



5^0 



RiBy 



Using the regulating function we can write Cj as: 



d=hme ijk L 3 [ J=!L=F* b {x) [ d 3 yf s (x,y)E al (y) [ d 3 zf 5 (x, z)E b \z). (15) 
JRi y/h(x) Jrj Jrj 



6^0 



To ensure proper contractions of the internal indices in the last expression we connect the points x, 
y, and z with holonomies of the Ashtekar connection along some smooth paths 71 and 72 connecting 
the points. We also use the identity e^k — — ATi[TiTjTk] to write Eq. ( |l5| ) as: 



2 The same result can be obtained || by imposing the gauge fixing at lagrangian level and after that 
performing the Legendre transformation. 
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Cj = - lim 4i 3 [ -JfL , 

■5^0 JRj VHx) Jrj 



(Tyf s (x,y) \ d 6 zfn{x,z) x 
Ri 



xTr[F Q A b (x)r fc C/ 7l (x, y)E m (y^U^iy, x)U l2 (x, z)E b > {z)t 3 U^ (z, x)]. 



7 2 



(16) 



To complete the regularization we also replace the curvature F^ b (x) by its approximation by a 
holonomy: 

1 



F« b (x)r k = lim— U(^ b] ), 



where 7 e - is a loop with area e 2 in the (a, b) coordinate plane, based at the point x. We have 
included explicitly the antisymmctrization with respect to a and b to ensure the vanishing of first 

x,[ab] 

d 3 x 



term in the expansion of the holonomy U{l c x r.g,) in powers of e 2 . Thus for Cj we get 



d d yf s (x,y) / d 6 zf s (x,z) x 



5"'5 2e 2 ./ H/ Jh{£)J Rl 



where we have used the convention 2_E al (y)r,; = E a {y). The expression under the trace is exactly a 
Smolin-Rovelli (l2| loop variable — T' Q& 1(7^ .jjz"/ xyx #-y xzx ) based on the loop shown on Figure [l]. 
Thus finally we get for the regulated version of Cj: 



where 




FIG. 1. Loop 7^ , a j,#7a;Hi#7a;za; on which the T-operator is based. 



Cj = lim lim C r 

e->0 5->0 ' 



d 3 :i 



d 3 yf 5 (x,y) d 3 zMx,z)Y,T [ab] (Y x&i #l X yx#J^)- (17) 
2e 2 J Rl y/h(x) J Rl Jrj 

ao 



Now we promote the last expression into an operator by replacing the loop variable with the 
corresponding loop operator. Thus for the Hamiltonian operator H we get: 



H = lim V . /- lim lim df Ae + AV? 



(18) 
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where Vi is the volume operator, as defined in pL As it was shown in |TT|] the spin network states 
are eigenstates of the volume operator so we can replace Vi with the corresponding eigenvalue for 
the volume of the I-th cube. 

Thus we have a regulated version of the Hamiltonian operator with which we act on the spin net- 
work states. We follow the standard procedure of regularization in which we apply the Hamiltonian 
on the states, perform all integrations, and at the end take the limits. 

Definitions of the spin network states can be found in H , || , [Q , . For our purposes it suffices 
to recall just the basic components of the definition. The spin networks are defined by: 

• a closed graph T in three-space; 

• labeling of the edges by irreducible representations of SU(2). We can interpret the labels as 
giving the number of loop segments along the corresponding edge; 

• intertwiners at the vertices, defining the way the loop segments coming from the edges are 
routed through the vertex. 

For a rigorous description of the way the spin networks are projected on a plane we would need 
some additional details but for simplicity in our calculations we will assume that the spin network 
we are using has been already projected. 

In our work we introduce a modification of the definition of the loop operators. In their standard 
definition the loop operators T ab [7] are based on a loop 7 and corresponding to every index there is 
a "hand" attached to the loop 7. In our case, because the loop on which T ab is based shrinks to a 
point , we have the freedom of modifying the attachment of the "hands" in a way convenient for our 
calculations. We consider the base loop 7 to be a planar loop with "hands" based on spin network 
edges of infinitesimal length i5 attached to 7. These edges have color 2 and are denoted j X yx and 
Ixzx in equation Eq. (|l7|). It will be also convenient for us to split the points at which the "hands" 
are attached to the loop 7 so there is a distance of order 5 between them. This can be thought of 
as a choice of decomposition of the fourvalent vertex positioned at the point x into two trivalent 
vertices. 

With such a definition it can be easily shown that the standard action of the "hands" is represented 
by connecting the spin network, on which T ah is based, by edges of color 2 to the original spin 
network. At the place of each grasping the action creates new trivalent vertices. Also we multiply 
by a factor of Wnlp^i A°[ei, 7(3)] whenever a "hand" situated at j(s), corresponding to a index "a" 
of the loop operator grasps an edge of color pi from the spin network. Ipi is the Planck length 
and A a [e,,7(s)] is the standard distributional expression: 

A a [e l;7 ( S )] = Idtei{t)6 3 (j(s),ei{t)) 
The action of one of the "hands" of the T a& -operator is shown on Figure ||. 




j,ab 

FIG. 2. Action of a "hand" of the f "^-operator. 



G 



Now we are in position to apply the Hamiltonian operator to the spin network states. We also 
assume that the spin network states are normalized as in [Tlfl . For the action of the Hamiltonian on 
(S, T\ we get: 

(S,T\Cf> 5 '< = £ [ J£= [ d 3 yf s (x iy ) [ d 3 zf s (x,z)x 
x E{E< e *l f[abI (<aS#7^#7*^ (19) 

ab i k i>J 

In the first term of the above sum i runs over all edges in the 7-th cube. In the second term the 
index k runs over all vertices inside the 7-th cube and the indices i and j run over all edges joined 
at the fc-th vertex. We get two types of terms: edge terms and vertex terms. In an edge term the 
two "hands" of the Hamiltonian grasp one and the same edge of the spin network, labeled e,, having 

(k) (k) 

color pi. In vertex terms, the two "hands" grasp two different edges, el and , joined at a vertex 
Vk ■ The two edges in general have different colors pi and pj . We divide the vertex terms further 
into two sub-cases: either the tangents to the edges at the common vertex are collinear, or there is 
some angle 9, different from 0° and 180° between the tangents. 

We will assume that the cubes have been shrunk enough so that in a single cube there is at most 
one vertex. The evaluation of the action in any case can be split into two different parts - we have an 
analytical part, coming from calculating the pre- factors in the action of the "hands" and evaluating 
the integrals, and a graphical part in which we complete the limiting procedure by shrinking the 
attached loops. 



III. ANALYTICAL ACTION OF THE HAMILTONIAN OPERATOR ON SPIN 

NETWORK STATES. 



A. Action on a single edge. 



For the first type of terms in the last bracket in Eq. (19) we get: 
L 3 f d 3 x 



2e 2 J Rl Jhfr 



f d 3 yf s (x,y) f d 3 zf s (x,z)Y,(e t \f lab] (^ ai> # 
Ri JRi ., 

ab 



lxyx#lxzx) = (20) 



— -fs{x,e % {s))fs{x,e l {t)) x 

y/h(x) 

J2 e t [a (s)e^ (*)(e i #^# 7xsB # 7x « c | 



x 



The meaning of the notation in (ei^f^ .•}H z ^xyxH z lxzx\ can t> e understood from Figure ^. The 
dashed circle denotes the region which will be shrunk to a point. To proceed further we have to 
specify the regulating function f$(x,y). We use normalized, weighted ^-function: 
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FIG. 3. Grasping of the Hamiltonian on a single edge. 



fs(x,y) = y / h(x)fs(x,y) = (-^ s )^/h(x)e[8 -\x-y\] 



(21) 



After evaluation of the space integral in Eq. (EG) we get 



Ri vHx) 



fs(x,ei(8))f s {x,ei(t)) 



8tt<5 3 



3d 
2~5 



(22) 



for d — d(s, t) — \el(s) — el(t)\ less than 25 and zero - otherwise. We have assumed also that 5 < L. 

As we are going to let 5 to go to zero, this will force also the separation between €i(s) and el(t) 
to approach zero. This is why we keep one of the parameters, say s, fixed and expand e'i(t) in 
power series about s. The first term in the expansion will make the whole expression vanishing, 
because of the antisymmetrization of the product e'J a (s)e'j b ' (s). Also the distance d can be replaced 
by \(t — s)\\ei(s)\. Combining Eq. (|2p| ) and Eq. ( p2|) we get to the lowest order in 5: 



L 3 3 
e 2 <5 3 167T 1 



(167rZ£j) 2 p? ds dt 



,|(i- S )||e;( S )| , \(t- 5)11^)1 



25 



25 



The limits of integration with respect to t are determined again by the expansion of d(s, t) and are 
given by: 



t e [s 



266- 



266+ . 



where 6 + — 1+0(6) and 6 — 1+0(8). As the integrals with respect to s and t are reparametrization 
invariant, we can choose a parametrization such that |e'j(s)| = 1. In this parametrization we set: 



e»(s) 



and 



e'i(s) 



P(s) 



where f and n are the unit tangent and normal vectors to the loop and p(s) is the curvature radius 
of the loop. Now we can perform the integration with respect to t and get: 
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r s+288+ 
ls-255- 

Thus for Eq. (§3) we obtain 



[-S+255+ 

/ dt\t - s\ 

Js-2SS- 



2-3 



26 \ 26 



6 2 {\ + 0{6)) 



Using the corresponding expression from connection representation in terms of holonomies it can be 
shown that up to terms of order e the last sum can be written as: 

^r[ a ( S )n b ]( S )(e 4 # 7 ^#7^ ;c #7- :c | = (ei#^ M #7*»*#7«*l (23) 

ab 

2 

where now the loop 7^ f^JFJxyxfFlxzx is in a plane defined by the tangent and the normal to the 
edge ej at the point s. Note that this loop is well defined as in the terms where the edge is a straight 
line and the normal is not defined, the curvature radius becomes infinite and such terms vanish. In 
the general case as the size L of the cubes goes to zero we can replace the integral with respect to s 
by its mean value: 

ds 2 Ki T j 2 

^y( e i#7l,f»#7^#7x Z x| = — (ei#7l,fn#7z ! /*#7x Z x|j ) 

where n is a number of order one, depending on the orientation of the edge inside the cube. We 
finally get the analytical expression for the action of the Hamiltonian on a single smooth edge: 

' ;V 3 ilMl^p^ltfn^y^lx^i (24) 



e 2 pi6 10ir 

This expression is a product of a pre- factor, containing a combination of regulating parameters and 
the original spin network with attached additional loops, subject to graphical evaluation. The pre- 
factor will determine the way in which we take the three limits e — ► 0, 6 — > 0, and L ~ * to make 
the whole expression finite. 

To understand the above intermediate result we compute the action of the regulated diffeomor- 
phism constraint: 

(7(iV) = hm VlimUm^E / d 3 xN a (x) [ d 3 yf s (x, y)^^*^] (25) 

j ab JHl JR ' 

on a single smooth edge of the spin network state. The result, up to finite numerical factors, is the 
same as Eq. (|24| ) if the shift vector N a (x) is given by: 

N a (s)=h a (s)-^-. (26) 
p{s)6 

This is not surprising since a similar result is true for the calculation in the connection rep- 
resentation. To obtain this result one also must perform the graphical evaluation of the loop 
7z fh^lxyx^lxzx from Eq. ( p4{ ) and of the loop 7^ fA #lxyx from the calculation of the action 
of the diffeomorphism constraint, which shrink at the end of the calculations. We show in detail this 
evaluation later in the paper in Section 4.1. 

Thus the action of the Hamiltonian on smooth edges of the spin network states can be interpreted 
as diffeomorphism transformation in a direction defined by Eq. (p6[). This means that if we act on 
diffeomorphism invariant states, the action of the Hamiltonian on edges annihilates them. As we 
assume that the states we act on a diffeomorphism invariant we will discard the action on smooth 
edges from the final result. 
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B. Action on two edges meeting at a vertex with collinear tangents. 




FIG. 4. Grasping of the Hamiltonian on two edges meeting at a vertex and having collinear tangents. 



An example of the grasping of the Hamiltonian on two edges meeting at a vertex and having 
collinear tangents is shown on Figure f|. The analytical calculation in this case is almost the same 
as with the case of a single edge. The only difference comes from the fact that the grasped edges 
could have different colors so instead of Eq. ( |24| ) we have: 

kL 4 3 2 
(16^) 

In this case we can not interpret the action of the Hamiltonian as a diffeomorphism transformation. 
As discussed in |l4| in such cases the action can change the diffeomorphism class of the spin network. 
But, as we will see later in the paper, such terms are of order L and they vanish after taking all 
limits. 



C. Action on two edges meeting at a vertex with non-collinear tangents. 



Let the Hamiltonian grasps the edges ej and ej with colors pi and pj respectively. Let and ej 
meet at some vertex Vk- We denote the part of the spin network state corresponding to the edges 
under consideration by ((ej, ej)|. For the action of the Hamiltonian we have: 



L 3 

((e t ,e 3 )\ — 



d 3 :i 



2e 2 Jr, ^Hx)J Ri 
L 3 



d 3 yf s (x,y) / d 3 zf s (i,z)f^ b \ 1 ^ ni # lxyx # lxzx ) 



2e 2 



{lQTTl I Pl Yp l Pj / ds / dt 



d 3 x 



--f5{x,e i (s))f 5 {x,e j {t)) x 



Ri \/H x ) 



ab 



Again ((ej, ej)#7^ & ^#lx V x4t z lxzx \ denotes the action of grasping of the Hamiltonian. Using Eq. ( p2] ) 
we can evaluate the spatial integral to get 

L 3 3 



e 2 S 3 16tt 



(167rZp;) piPj ds dt 



2-3 



|e1(s) - ej(t)\ f\ei(s)-ej(t)\ 



26 



2<) 



(28) 



ab 
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We use a parametrization such that \el\ = |ej| = 1 so the distance d(s, t) = \el(s) — e}(t)| becomes 
simply d(s, t) = \s — 1\. Also let r] and C be the unit tangent vectors at the vertex and let the angle 
between them be 0y . Performing the integration with respect to s and t we get for Eq. (^8|) : 

1*8 loS^ 167 ^ 2pift ? ^ <(ei ' ^)#<^#7^#7^|/. 

Again we can use a formula which can be proven by using the corresponding expression in the 
connection representation, namely: 

J2 V la C b] ((ei, e J -)#7^ ag #7^#7^|i = sm0((ei, e^#7^#7*s/*#7^|i, 

a 6 

to transform the sum in Eq. ( p8| ) . In the last expression 7^ - is a loop in the plane defined by the 
two tangent vectors 77 and Thus finally we get for the action of the Hamiltonian on a vertex: 

L 3 38, 



e2§ 107r (16tt/ k ) PiPj{{euej)#7l A ^ X y X #7 xzx \i. (29) 

In the obtained expression we can separate between the pre-factor and the loop-deformed spin 
network state. The pre-factor contains a combination of the regulating parameters and an implicit 
dependence on the arbitrary background metric through the angle 6ij. The loop deformation will 
be subject to a graphical evaluation later in the paper. 



D. Taking the limits. 

Now we are in position to take the limits in the computation. Basically the limits appear in our 
calculation in two different ways. First the pre-factors in Eq. ( p7| ) and Eq. (^9|) contain combinations 
of the regulating parameters. The limits should be taken in such a way so these combinations yield 
finite results. Second, the loop on which the loop operator T ab is based shrinks to a point, together 
with its "hands" and this leads to the graphical evaluation. 

The combinations of parameters are: 

kJJ^ L 3 mr\\ 

^ and (30) 

for Eq. ( p7|) and Eq. (|2^ ) respectively. To insure finiteness of our expressions let us take the three 
limits e — > 0, S — > 0, and L — > along a plane in the (L, e, S) parameter space, defined by: 

where Z is an arbitrary constant, chosen in such a way so the relative order of taking the limits is 
satisfied. This order is determined by the conditions that 5 < L and 5 <C e. We have to be careful 



about the way the limits in L and 5 are taken, as the diffeomorphism vector Eq. (26) contains the 
ratio L 3 /S, so we would not like 6 to go to zero much faster than L. If we set LrjE — const then 
this will be enough to prevent the diffeomorphism from becoming infinite. 

We can now write the general formula for the analytical action of the Hamiltonian on diffeomor- 
phism invariant spin network states. All the cubes which are empty give zero. 

The result from the action on different edges meeting at a vertex with collinear tangents Eq. ( |27| ) 
can be written as : 



3kZL 
lOnpi 
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and vanishes as L goes to zero. 

Adding up all contributions from the vertices with non-collinear tangents we get: 



(S,T\Cf^ = _(l 67 r/|,0 2 ^^%((ef\ef ) )# 7 ^#7^#7 :cz:c | (31) 

i,3 

In the above sums the indices i and j run over all edges joined at the vertex which is in the 7-th 
cube. The sum contains the arbitrary but finite constant Z and the angle 6ij between the edges at 
the vertex Vk ■ The bar over the state represents the fact that we still have to perform the shrinking 
of the attached loops and thus to evaluate the graph. 

We get overall action which is finite but background dependent because of the presence of the 
angle 6ij . This is a problem which can be solved by redefining the way we approximate the curvature 
F a i,. Let the holonomy which approximates the curvature be based on a loop with area Oije 2 
instead of e 2 . Then all the calculations go through but at the end there is an extra factor of 
Oij in the denominator of Eq. ([3l| ) to cancel the corresponding factor from the numerator. It is 
important to notice that although the explicit inclusion of the factor of 6^ formally solves the 
problem of background dependence, the situation is not completely satisfactory. The computation 
of the evolution of the spin network states involves at each step measuring the angle Oij between 
each pair of tangent vectors. Technically this might require the introduction in the definition of the 
Hamiltonian of an operator, measuring the angle 6^ . This issue requires further investigation. 



IV. GRAPHICAL ACTION OF THE HAMILTONIAN. 

To evaluate the deformed spin network state from equation Eq. ( pT| ) we use some techniques from 
the recoupling theory of colored trivalent links and knots. First let us consider a general situation 
in which a loop attached to a spin network is shrunk to a point. An important observation is that 
the evaluation is local, in the sense that whatever result we obtain, it is independent of the way the 
spin network is connected outside of the circle denoting the shrinking region. 

To evaluate the attached loop we have to do the following 

• Expand all the edges which are entirely in the shrinking region as an antisymmetrized sum 
with all possible crossings between the loop segments. 

• Resolve each crossing of single loop segments in the shrinking region according to the binor 
Mandelstam identity (see Q for discussion). 

• Associate to each closed single loop the "loop value" (-2). 

• Smoothen the obtained graph to get if possible the original spin network. 

But these are exactly the operations which occur in the evaluation of the Kauffman bracket |HJ 
. The Kauffman bracket is an invariant of regular isotopy of colored knots and links with trivalent 
vertices. Thus for the evaluation of the graphical action of the loop operators we can use the 
techniques of recoupling theory, developed for computing the Kauffman bracket pl|JT^ -p7[ . What 
makes these techniques powerful is the fact that we do the computations blockwise. First we work 
out the analytical expressions for the most simple graphs. Then in more complicated calculations 
we identify the simple graphs and replace them with the corresponding analytical expressions. This 
is possible because the evaluation of the graphs is local - we can shrink to a point just a portion 
of the whole graph, keeping the remaining parts fixed. The basic formulae of the recoupling theory 
are summarized in the Appendix. With the use of these formulae we compute the graphical action 
of the Hamiltonian on smooth edges and on vertices. 
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A. Evaluation of the graphical action for edges. 



We perform the graphical evaluation for the grasp of a smooth edge to show that the graphical 
action of the Hamiltonian on smooth edges is equivalent to diffeomorphism transformation. We start 
from the expression (ei#j x f fl 4t z lxyx z tf z lxzx\- The grasping is shown on Figure ||. First we use the 
basic recoupling formula Eq. (Al) from the Appendix on one of the edges. As a result a new link 
of color n appears in the new configuration. The allowed values for n are determined by the basic 
properties of the trivalent vertices - the only values it can take are pi ± 1. 



(ei#lx,fn#Jxyx#Jx 




El Pi Pi n 
1112 

n— pizkl 




On the right hand side appears the 6j-symbol, as defined in |U|. We repeat the same step with the 
other "hand" and get: 



n—pi ± 1 m— pi±l 





Pi 


n 1 


/ Pt 


Pi 


777 1 




1 


2 J 


X 1 


1 


2 J 



Pi! n Pi m \ Pi 



As a next step we remove the internal "bubble" using the i?-net from equation Eq. (A4) from the 
Appendix to obtain: 



XZX | 



\- U n f Pi Pi n\ #(pj,l,n) EIL 
' I 1 1 2 / n + l 




To continue further we have to make more careful analysis of the way the limit is taken, when the 
states are diffeomorphism invariant. As this is not relevant for our considerations we will stop here 
and show that the above expression is analogous to the result from the action of the diffeomorphism 
constraint. To this aim we act with Eq. ( p5| ) on a smooth edge. There the loop deformation is given 
by the expression 7^ ffi ^ xyx . To evaluate graphically this expression we use the basic recoupling 



formula Eq. (Al) and immediately get: 



{ e i4^1x,fh 



i^^lxyx I 




El Pi Pi n 
1112 

n— Pi±l 




Thus we see that in both cases of the action of the Hamiltonian and of the diffeomorphism con- 
straints, the graphical evaluation leads to the a "bubble" on the smooth edge with different numerical 
factors. This proves that the two graphical actions are equivalent and thus we can disregard all terms 
in the action of the Hamiltonian on smooth edges of diffeomorphism invariant states. 
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B. Evaluation of the graphical action for vertices. 



We apply the techniques from the recoupling theory of colored graphs also to evaluate the ac- 
tion of the Hamiltonian when it grasps two edges meeting at a vertex. For simplicity we consider 
only the action on one pair of edges joint at a trivalent vertex. We have initially the expression 
((ei,e :) )#7 e .^jxyx^Jxzxl which can be transformed with the use of the recoupling formula Eq. 



(Al), given in the Appendix. We obtain the following result: 




El Pi Pi n 
1112 

n— pi±l 




Then we repeat the same step with the other "hand" of the T-operator to get: 



E E 

n— pi ±1 m—pjzLl 



Pi Pi n 
1 1 2 



Pi Pj m 
1 1 2 




The the inner triangular diagram of the above graph can be evaluated with the use of one recou- 
pling formula and one "bubble" removal (see the Appendix): 

((e i)e ,)#7^#7^#7^l= E E 



n p t p k \ (-l-) Pk &(pi,Pj,Pk) 



Pj m 1 



Pk + 1 




The last step in the graphical calculation requires more careful consideration. As discussed in || 
when we act on diffeomorphism invariant states, the spin network remains in its knot class, although 
the area of the last remaining loop goes to zero. To understand this situation let us discuss the action 
of the Hamiltonian from the point of view of diffeomorphism invariance. We start with a state which 
is diffeomorphism class of spin network and a Hamiltonian, also based on a class of diffeomorphism 
invariant smooth loops with "hands" . To perform the action we introduce a background metric, 
break the diffeomorphism invariance, and introduce notions of lengths and areas. It is in this non- 
invariant sense in which we can talk about "loops of some area shrinking to a fixed point" and about 
"hands of infinitesimal length" . Also all the formulae we apply from the recoupling theory are based 
on the ideas of replacing loops by their loop values and recoupling, which are diffeomorphically 
non-invariant operations. We assume that we can perform these operations as far as we recover at 
the end the "right" diffeomorphism class. What we expect at the end is the original spin network 
attached through two "hands" to a loop which in the non-diffcomorphism limit shrinks to a vertex. 

We can choose the attached loop in different way but we want our choice to be consistent with 
the operations we perform in a non-invariant fashion. One possible such choice is the attached loop 
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to connect two edges joined at a vertex and then to run along the edges of the original spin network. 
This means that when the state we are acting on is diffeomorphism invariant taking the last limit is 
trivial - the graph does not change. Thus finally we get: 



(S,T\Cj' 5 '' 



3Z 
IOtt 1 



(16tt1 2 pi ) 2 2_ j p l p j 



i,3 





Pi 


m 1 


( n 


Pi 


Pk \ 


{7 


1 


2 J 


I Pi 


m 


1 J 




(32) 



(-l) p *tf(p«,Pj,Pfc) 
Pk + 1 



where the vertex having the additinal edge attached to it is in the 7-th cube. The state (S, T\ we 
acted on is normalized according to the normalization introduced in p"T[ . If the final states are also 
to be normalized we have to introduce some additional factors in Eq. (|3^). As the old trivalcnt 
vertex has been transformed into a new one, two new edges and two new vertices have been added 
we get for the final normalized sum of states: 



(S,T\Cf 



Pi Pi n 
1 1 2 



LAe _ 3Z(16irl 2 pl ) 



IOtt 



Pj Pj m 
1 1 2 



E 



PM z2 E 



I -2A n A m i9(pi,pj,pk) 
#(l,Pi,n)#(l,Pj-,m)i?(7i,m,pfc) 



(33) 



n pi p k 
Pj m 1 



-l) p *tf(pi,p 3 -,Pfc) 
Pk + 1 




The action of the Hamiltonian constraint C as part of the physical Hamiltonian H can be described 
as follows: When the Hamiltonian constraint C acts on edges, the action can be interpreted as a 
diffeomorphism transformation. Thus the Hamiltonian constraint acts non-trivially on diffeomor- 
phism invariant states only when it acts on the vertices of the spin network. The action on vertices 
amounts to adding a new edge attached through new trivalent vertices to each pair of edges at each 
vertex. Each new edge has color 1. Also the colors of the edges connecting the new vertices and the 
original one change. This change can be explained with the requirement that the new graph is again 
a spin network. The obtained states are multiplied by finite factors which, although cumbersome, 
are straightforward to compute. 



V. THE HAMILTONIAN EIGENSPECTRUM 

Let us recall that the Hamiltonian defined through the "clock" field has the form Eq. (fl8|): 



H = lim V J- lim lim df' 5 ' e + AV?. (34) 

To find the eigenspectrum of the Hamiltonian operator we have to handle in a satisfactory manner 
the square root in Eq. (|34|). Since we have not been able to complete this task we will only describe 
the directions for future work. 

Since both the Hamiltonian constraint and the volume operator in Eq. ( |34"| ) give non-zero results 
only when acting on vertices, the above sum reduces to a sum only over the cubes in which there is 
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a vertex. Thus the sum in Eq. (|3J) becomes finite. Also all the terms coming from the action on 
different vertices are independent from each other. This justifies the following strategy for computing 
the Hamiltonian eigenspectrum: We can use the result from the action of the expression under the 
square root on vertices - the action of the first term is given by Eq. (Q) and the action of the volume 
piece have been computed for example in pLlj ] . To proceed further we might have to distinguish 
between the action on bivalent and trivalent vertices on one hand, as having zero volume, and the 
action on higher- valence vertices. In any case we will have to diagonalize the expression Eq. (|33|). 
The problem is that a priori we can not expect the matrix defined by Eq. ( |33|) to be diagonalizable. 
Since the expression for the Hamiltonian constraint we start with is a real one we can assume that 
there exists an appropriate symmetrization which can make the constraint operator self-adjoint and 
its matrx symmetric. This is why to make Eq. ( |33| ) diagonalizable we could replace it by the semi- 
sum of it and its transposed. Then we will have to diagonalize the sum of the Hamiltonian and 
volume matrices and to take the square root. The result will be an infinite dimmensional diagonal 
matrix with the Hamiltonian eigenvalues as its elements. 



VI. CONCLUSION 



The action of the Hamiltonian constraint on a spin network state can be described as a transition 
of the spin networks on which the state is based from one knot class of graphs with vertices into 
another one. In this transition the coloring of certain edges changes so the requirements for the 
graph to be a spin network remain satisfied. At the same time each new state is multiplied by a 
factor, which carries the appropriate dimensions, contains an arbitrary constant, and also includes 
the result from the computation in recoupling theory. 

In the paper we considered in detail only trivalent vertices but the generalization to higher valence 
vertices is straightforward. As discussed in |fjc[ , [jy], any higher- valence vertex can be decomposed 
into a set of infinitesimally displaced trivalent ones. The analytical calculation we performed did 
not depend to the valence of the vertex. The difference with the considered case will appear in the 
graphical evaluation where we again can apply the recoupling theory. It can be shown that as a 
result in the same way as with the trivalent vertices, edges of color one are added, re-routing of the 
loops through the vertices occurs, and one can compute the corresponding pre-factors. 

Also the evaluation presented in the paper can be applied to the q-deformed spin networks Jl6| |, 
. In the q-deformed theory all the formulae we used have their counterparts. Thus in that case 
one should simply replace the formulae we derived with the corresponding q-deformed versions JL5| . 

As we discussed before, to get rid of the background dependence of the pre-factor we have to know 
the angle between each pair of tangents, which requires modification of the regularization. Another 
set of problems occurs when we try to make sense of the square root involved in the definition of the 
Hamiltonian. There are certain proposals in this direction: One |l§| ] is to use an expansion in inverse 
powers of A and to find an approximation to the square root. Another proposal j(| involves different 
procedure, which also allows the square root to be expressed as a series. In any case we have to 
settle this set of problems before trying to solve for the eigenvalues of the Hamiltonian operator. 



This work would be impossible without the numerous discussions I had with Lee Smolin and 
Carlo Rovelli. I also would like also to thank Don Neville and Seth Major for their comments and 
criticism and the Center for Gravitational Physics and Geometry where this work was started, for 
the hospitality. 
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APPENDIX: BASIC FORMULAE OF THE RECOUPLING THEORY. 



One of the main results of the recoupling theory of colored knots and links with trivalent vertices 
[fL5[ is the computation of the Kauffman bracket for different framed spin networks. The framing 
refers to the fact that in the computations one keeps track of over- and under-crossing. In our work 
we do not make this distinction so we use a simplified version of the recoupling theory, namely we 
replace the deformation parameter q by its "classical" value —1, relevant to our case. As a result 
in most of the formulae from [ fL5| the q-deformed (or quantum) integers are replaced by ordinary 
ones. We list here the basic formulae, which we use, with the mentioned corrections made. For more 
details see p5| . 

• The basic relation in this theory expresses the relation between the different ways in which 
three angular momenta, say j\ , and j% can couple to form a fourth one, j'4. The two possible 
recouplings are related by the formula: 




where on the right hand side is the g-6j symbol, as defined in |15(| for the value of the defor- 
mation parameter q equal to -1. Again the dashed line denotes the fact that the recoupling 
occurs in a region which shrinks to a point; it is not extended in space. 

• Closed loops which have been shrunk to a point are replaced by their loop value, which is (for 
a single loop with zero-self-linking) equal to —2. The evaluation of a single unknotted loop 
with color n is fllal: 



(-1)> + 1). (A2) 



The small rectangle in the above diagram denotes the antisymmctrization of the n-line. 

• The next graph we consider is the "bubble" diagram. Upon shrinking the "bubble" , this 
diagram will reduce to a single edge so the evaluation will be different from zero only if 
the colors of both ends of the "bubble" are the same. Thus the "bubble" diagram equals 
some numerical factor times a single edge. By closing the free ends of the diagram it is 
straightforward to show that: 



o 

b 



= Sr, 



(-l)"0(a,6,n) 
n + 1 



(A3) 



in which the function 6(a, b, n) is given, in general, by 
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6{m, n, I) 




= , n (a+6+c) (a + b + c+iy.a\b\c\ 
1 ' {a + by.(b + c)l(a + c)\ 



(A4) 



where a + 6 = m, a + c = n, and b + c = I. 

A basic element in most of the formulae of the recoupling theory of colored graphs is the 6j 
symbol. It is defined through the so called Tetrahedral net via the following relation 



a b e 
c d f 



(-l) e (e + l)Tet 



a b e 
c d f 



i?(a, d, e)i?(c, 6, e) 



(A5) 



The Tetrahedral net is represented by the following diagram 




Tet 



a b e 
c d f 



Upon evaluation the Tetrahedral net yields: 

1 



Tct 



a b e 
c d f 



E 



£ ^ TT (s-ai)! 11 (bi-s)\ ' 



where 



ai = ±(a + d + e), &i = \{b + d+e + f), 

a 2 = ^(fe + c + e), b 2 = |(a + c + e + /), 

a 3 = |(a + &+/), 6 3 = i(a + 6 + c + d), 

a 4 = i(c+d + /), 

m = max{ai}, M = min{6j}, 

£ = a! 6! c! d\ e! /!, X = Unfa - <n)\ . 



These formulae are sufficient for the computations performed in the paper. In |lq| one can 
find a detailed derivation of all of them. 
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